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ABSTRACT 



We study the distribution of stars, HII regions, molecular gas, and individual 

w ■ giant molecular clouds in M33 over a wide range of spatial scales. The clustering 

strength of these components is systematically estimated through the fractal 

dimension. We find scale-free behavior at small spatial scales and a transition 

to a larger correlation dimension (consistent with a nearly uniform distribution) 

at larger scales. The transition region lies in the range ~ 500 — 1000 pc. This 

transition defines a characteristic size that separates the regime of small-scale 

turbulent motion from that of large-scale galactic dynamics. At small spatial 

scales, bright young stars and molecular gas are distributed with nearly the same 

three-dimensional fractal dimension (-D/,3D fS 1-9) 5 whereas fainter stars and HII 

£T) ■ regions exhibit higher values -D/,3D — 2.2 — 2.5. Our results indicate that the 

£->! . interstellar medium in M33 is on average more fragmented and irregular than in 



the Milky Way. 

Subject headings: galaxies: individual: M33 - - galaxies: structure - - stars: 
formation 



1. Introduction 

In the Milky Way, gas and dust are organized in a hierarchi cal and self-similar manne r 



that it is supposed to be a consequence of turbulent processes (jElmegreen fc Scald 12004 ). 
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These fractal patterns are observed over at lea st the range 0.1 < r 5, 100 pc, going from 
dense cores to giant molecular clouds (GMCs, iBergin &: Tafallal 120071 ) . The formation of 
stars also shows fractal features, usually observed in regions with a spatial hierarchy ranging 



from a few pc up to about a kpc for so-called star complexes (e.g. lEfremovlll995l ; lElmegreen 



2010l ). Star clusters are in the lower levels of this hie rarchy, although there is evidence that 



youn g open clusters also exhibit smaller substructure (ISchmeja et al. 
2009h . 



2008; Sanchez fc Alfaro 



This complexity can be characterized in many different ways. Some common strategies, 
such as estimating mass or size distributions of certain types of objects (clouds, cores, clus- 
ters), have to be taken with extreme cauti on because they depend, among other thi ngs, on 
the criteria adopted to define the objects (IPineda et al.ll2009l ; ICurtis &: Richer! 120101 ) . Note 
that in a rigorously hierarchical scenario there is no characteristic spatial scale that can be 
used to define any particular structure. Moreover, it has been shown that pr ojection effects 



can significantly alter t he estimated masses and sizes of molecular clouds (ISanchez et al. 



20061 ; IShetty et al.ll2010l ). Many other tools are now wi dely used to describe the complexity 
of these structures objectively and quantitatively (see lElmegreen fc Scald 12004 for a com- 
prehensive review). Fractal analysis is particularly appropriate for dealing with hierarchical 
and self-similar systems. The fractal dimension Df, which quantifies the spatial heterogene- 
ity, can be calculated for both the distribution of gas and the distribution of star-forming 
sites. This approach allows for the comparison of similar objects in different regions and/or 
under different physical conditions, and also the comparison of different types of objects. 

It is often accepted that the fractal dimension observed i n the interstellar medium (ISM ) 
has a nearly universal value around Df ~ 2.3 ± 0.3 (e.g. lElmegreen fc Falgarond Il996l ). 
This universality would indicate either that interstellar turbulence is driven by the same 
physical mechanisms everywhere, or that different physical mechanisms can (and tend to) 
generate essentially the same type of structures. However, the robustness of this conclusion 
is questionable given the wide variety of results reported and their associated uncertainties. 
If Df is inferred from prop erties such as cloud masses or sizes, the resulting uncertainties 
may be unacceptably large (ISanchez et al.ll2006l ). A more appropriate strategy is to measure 
the fractal dimension directly from an observed map. The boundaries of the projected 
images of interstellar clou ds in the Galaxy ha ve fractal dimension values spread over the 
range 1.2 < D pcr < 1.5 (ISanchez et al.l 120051 ) but it is not clear whether there are real 
variations from region to region or whether the different values reflect different observational 
data and/or analysis techniques. Moreover, even though D pcr ~ 1.3 is a constant, the 
assumption that one can relate D pcr to the fractal dimension in three dimensions using 



Df = Dpcr + 1 — 2.3 has been shown not to be valid (ISanchez et al.ll2005l ) 
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In a previous work we analyzed several emission maps of three different molecular clouds 
(Ophiuchu s, Perseus and Orion) and obtained Df ~ 2.6 ±0.1 with no evidence of significant 
variations (JSanchez et al.ll2007bl ). Similar fractal patterns should be observable in the distri- 
bution of newborn stars if the distribution of high-density cores follows the spatial structure of 
the parental cloud. In fact, the young mas sive stars in the Gould Belt exhibit a fractal pattern 
with the similar value Df = 2.68 ± 0.04 (jSanchez et al.ll2007al ). However, a recent study by 
Schneider et al.l (120101 ) reports remarkable differences in the A- variance spectra between low- 
mass star forming clouds and massive GMCs. Star clusters also show a wide variety of spatial 
patterns even for young, e mbedded clusters (ILada k Ladal 120031 ; ICartwright k Whitworth 
2004 ; ISchmeja et al.ll2008l ). but it is not clear whether thes e variations are due to evolu- 



tion or to differences in t he structure of the original clouds (JGoodwin k Whitworth! 12004 
Schmeia k Klessenl liooS I Allison et al"1l2009h . 



Simulations of turbulent fluids produce very diffe rent structures d epend ing on which 
processes are considered in the system. For example, iFederrath et al.l ( 120091 ) showed that 
simulations of supersonic isothermal turbulence in the extreme case of purely compressive 
energy injection, produce a significantly smaller fractal dimension for the density distribu- 
tion (Df ~ 2.3) than in the case of purely solenoidal forcing (Df ~ 2.6). Although it is 
widely accepted that turbulence is the primary driver of the structure and motion of the 
ISM, the main energy sources for this turbulence are not yet well established. Part of the 
problem lies in the wide variety of physical processes that can generate turbulent motions, 
such as protostellar jets and outflows, (proto) stellar winds and ionizing radiation, expanding 
HII regions, supernovae, cloud collisions, galaxy interactions, and gravitational, magnetoro- 
tational and other fluid instabilities (recent reviews on interstellar turbulence can be seen i n 
Elmegreen k ScalclliooJ : iMac Low k Klessenl 12004 : iBurkertJliooa llVfcKee k Ostrikerl 120071 ) . 
It is reasonable to expect that, depending on the dominant physical mechanisms driving the 
turbulence, the resulting structure may differ from region to region in the Galaxy. 

A related issue is the spatial extent of this self-similar behavior, since different physical 
mechanisms might dominate on different scales. In the solar neighborhood, fractal behav- 
ior has been observed fo r the distribution of young open cluster and young s tars at spatial 
scales of up to ~ 1 kpc (Ide La Fuente Marcos k de La Fuente Marcos! 120061 : ISanchez et al. 



2007al : Ide la Fuente Marcos k de la Fuente Marcos! 120091 ). In external galaxies, hierarchical 



structures extend up to > 1 kpc scales for the gas and for stars and star- forming sites (some 



recen t examples are in iBastian et al.l 120071: iDutt a et al.l l2008t iGieles et al.l l2008t lOdekon 



2008 



2009 



Sanchez k Alfaro 12008: Bastian et al. 



2009 



Dutta et al.l l2009bl : IScheepmaker et al. 



Bonatto k Bicall2010l ). Som e galaxies show a change in clustering strength with scale. 



This has been seen in both g as (JElmegreen et al 



20071 : iDuttaet al.ll2008l . l2009al ) and young stars (lOdekon 



2001 



Padoan et all 1200 ll : iKim k Park 



2008). There is also evidence for 
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variations among galaxies. Recent work indicates that the power spectrum of the gas distri- 



bution is steeper for gal axies with a larger star formation rate per unit area (jWillett et al. 



2005c iDutta et al.ll2009bl ) . This trend is consistent with the claim that bright galaxies tend to 
distribute their star-forming s ites in less clustered patterns (IParodi fc Binggelill2003l ; lOdekon 
20061 ; ISanchez fc Alfaroll2008l ). In other words, a larger star formation rate in a galaxy tends 
to be correlated with a larger fractal dimension. Thus, at least on galactic spatial scales, a 
universal picture for the fractal properties does not seem to hold. 

A significant challenge in interpreting different measurements is that authors present 
fractal dimensions for different ranges of scales, identify their samples in different ways, and 
make different assumptions regarding the comparison between two-dimensional measures of 
clustering strength and the implications for the actual three-dimensional structure. Our 
approach here is to consider a case study in which we systematically analyze the clustering 
of different components of a single galaxy over a wide range of spatial scales. Because of its 



proxim ity, large size, and low inclination, M33 is a suitable object for this task. iBastian et al. 



( 120071 ) showed that the star formation in this galaxy appear s to be i n heren tly hierarchical, 
with no characteristic size for star-forming regions. However, lOdekonl (120081 ) found a transi- 
tion on a scale of several hundred pc in the slope of the autocorrelation function for young 
stars. In this work, we study the clustering strength in the distribution of young stars, 
HII regions, molecular gas, and individual giant molecular clouds. Treating each of these 
components as consistently as possible, we can directly compare the clustering strength of 
each component as a function of scale. In addition, we examine in detail the steepness of 
the transition in clustering strength, evaluate the range of scales over which it occurs, and 
analyze its meaning in the context of the maximum spatial scale of coherent star formation. 
We apply robust algorithms that have been previously developed and tested on real data and 
simulated fractals. Section [2] describes the data and Section [3] describes our calculation of 
the fractal dimensions. An analysis and discussion of the results are presented in Section |U 
Finally, a summary and the main conclusions are given in Section 



2. The data 



The catalog of lMassey et al.l (120061 ) . available via the VizieR[j database, provides UBVRI 
photometry of 146, 622 stars in M33 down to apparent magnitudes of 23 with photometric 
erro rs less than 10% . From this catalog we use the same set of young stars (ages < 30 Myr) 
that lOdekonl (120081 ) used to calculate the angular two-point autocorrelation function, that is, 



1 http://vizier.u-strasbg.fr 
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stars fulfilling the criteria —0.3 < V — I < 0.0 and —6.0 < Mj < —4.0. This set of stars was 
divided into two subsets that we refer to simply as "bright" stars (—0.3 < V — I < —0.1 and 
-6.0 < Mr < -5.0) and "faint" stars (-0.3 < V - / < 0.0 and -4.5 < M 7 < -4.0). The 
total numbers of stars for the bri ght and faint sets a re 534 and 1644, respectively. For the 
HII regions we use the catalog of iHodge et al.l (119991 ) . which gives the positions of 1272 HII 
regions in M33. From this total data set, we remove regions classified as unresolved, diffuse, 
linear and/or any other factor that may raise doubts on the real nature of such regions. We 
also remove the regions having null integrated Ha fluxes in the catalog. This means that we 
only consider regions having an isophote level of at least 10~ 16 erg s _1 cm -2 arcsec" 2 . With 
these requirements the total number of "bright" HII regions to consider i s 617. The distri- 



bution of giant molecular clouds in M33 is obtained from the catalog of iRosolowsky et al. 



( 120071 ) (also available through VizieR) from which we extract the positions of 149 GMCs. 



We adopt a position angle of 23 degrees and an inclination of 55 degrees (taken from 
Hyperledao) to deproject the positions of stars, HII regions and GMCs. This is a first step 
necessary to avoid the calculated fractal dimension becoming smaller than its true value 
( JSanchez fc Arfaroll2008). To convert the angula r sizes into linear sizes, we assume a distance 
of 960 kpc ( JBonanos et al.ll2006t lU et al.l 120091 ) . The positions of bright stars, HII regions 
and GMCs relative to the galactic center aj2ooo.o = l h 33 m 51 s , <5j2ooo.o — +30°39'37") are 
shown in Figure [TJ 

As an additional way to study the the distribution of molecular gas, we also directly use 
the high resolution maps of C O (J = 1 — > 0) emission for the center region of M33, kindly 
provided by Erik Rosolowsky ( IRosolowsky et al.l 12007). To avoid problems that might arise 
when estimating the fractal properties in noisy maps (jSanchez et al.ll2007bl ). we use only the 
combined NRO+BIMA+FCRAO data cube which has the highest signal-t o-noise ratio. The 



noise varies across the map but the typical rms noise temperature is 60 mK ( IRosolowsky et al. 



20071 ). The final resolution is 20" (93 pc). We collapse the data cube to produce the map of 
integrated intensity of CO emission shown in Figure [2j 



3. Estimation of the fractal dimension 



The degree of clustering of a point distribution can be quantified by the correlation 
dimension D c . For scale-free sets, this can be calculated from C(r) ~ r Dc , where the cor- 
relation integral C(r) is the average number of stars within a dista nce r of e ach st ar. In 
practice, the distribution of stars is not truly scale-free. For example, lOdekonl (120081 ) found 



2 http://leda. univ-lyonl.fr 
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Fig. 1. — Spatial distributions of (a) bright stars, (b) bright HII regions and (c) GMCs in M33. 
The axis coordinates are positions relative to the galactic center in kpc. The inset in panel (c) 
shows the area corresponding to the molecular gas map (Fig. [2]). 



a transition in the clustering strength of young stars in M33 at a scale of about 300 pc. 
Because the correlation integral for any particular scale r includes information about the 
clustering on all scales less than r, any deviation from a scale-free distribution will affect, 
to some degree, the form of C(r) on scales larger than the deviation. For this reason, it is 
sometimes preferable to use the differential form c(r), where C(r) ~ J c(r')dr', especially 
if the goal is to determine the scale at which a change in clustering strength takes place. 
Because of its differential nature, however, c(r) is too messy for small data sets. 

Another reason the distribution of stars is not actually scale-free is the finite size of a 
galaxy. For an otherwise scale-free distribution with a well-defined edge, one approach to 
dealin g with this lim it ation i s to include, for each scale r, only those stars at least r from the 
edge. ISanchez et al.l (j2007al ) used the minimum-area convex polygon that contains all the 
points of finite two-dimensional fractal distributions, and showed that it is an effective way 
of performing an edge correction. We use this method to create an edge-corrected version for 
C(r). Finally, on small scales the observed distribution of stars is not scale-f ree because of 



resolut ion limits and, ultimately, because of the size of the stars themselves. ISanchez et al. 



( J2007af ) found that the calculation of C(r) for two-dimensional fractals is reliable on scales 
for which the standard deviation of C(r) is not greater than C{r) itself. We use this as our 
guide in determining the smallest scales to include in our fits to the data. 



Figure [3] shows our results for correlation integral C(r) and the differential form c(r) 
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Fig. 2. — Integrated intensity of CO emission in the central region of M33 (Fig. [T]). The gray scale 
runs linearly from (white) to the maximum value (black, ~ 7.8 K km s _1 ). Contour levels are 
shown at 1 and 2 K Km s _1 . 



for each type of point distribution in our sample: bright stars, faint stars, HII regions, 
and GMCs. Each panel includes three functions; from top to bottom, they are the edge- 
corrected correlation integral C(r)/r 2 , the uncorrected version of C(r)/r 2 , and the differential 
correlation function c(r)/r. The functions are normalized so that a random two-dimensional 
distribution (corresponding to D c = 2) would produce a horizontal line. In other words, the 
slope on the log-log plot is D c — 2, so that a more strongly clustered distribution decreases 
more quickly with scale. A comparison of the bottom two functions in each panel illustrates 
the fact that the differential form c(r) is more sensitive to changes with scale, but messier, 
than the integral form C(r). A comparison of the top two function in each panel illustrates 
the effect of the edge-correction algorithm, flattening the slope on large scales. 

Gray shade d line segments in Fig ure p show our power law fi ts to the edge-corrected 
C(r). Following ISanchez et al.l (j2007al ) and I Sanchez fc Alfarol (120081 ). the range in r used for 
the fits is limited on small scales to regions where the variation in C(r) is not greater than 
C(r) itself; these lower limits are ~ 100 pc for stars and HII regions and ~ 450 pc for GMCs. 
In order to consider the possibility of a transition in D c , we performed separate fits at large 
and small scales. The range of spatial scales for the transition was varied and the result that 
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Fig. 3. — From top to bottom in each panel, the edge-corrected correlation integral C(r)/r 2 , the raw 
correlation integral C(r)/r 2 , and the differential correlation function c{r)/r for the distributions of 
pointlike objects. The data have been arbitrarily shifted vertically for clarity. Shaded line segments 
show power-law fits. 



minimizes the sum of the squared residuals, is that shown in Figure [3] and listed in Table [TJ 
We calculated uncertainties in D c , also listed in Table [JJ using a bootstrap algorithm. 

All of the objects except the GMCs exhibit scale-free clustering on small scales and a 
clear transition to a natter slope near 1 kpc. The transition region is found to lie in the 
range 500 < r < 1000 pc for bright stars and 600 < r < 950 for faint stars and HII regions. 

We use a very different type of clustering measure on the projected CO emission map 
(Fig. [2]). For th ese data, we ca l culate d the perimeter-based dimension D per using the same 
algorithm as in ISanchez et al.l (120051 ). We chose brightness levels from the lowest to the 
highest values in steps of 0.5 K Km s _1 . At each brightness level the algorithm defines 
"objects" in the image as sets of connected pixels whose brightness value is above this 
level. Then, the perimeters P and areas A of each object are determined and the perimeter 
dimension is obtained from the relation P ~ ^^W/ 2 , yy e ^j^ no ^ consider objects containing 
l ess than 20 pixels because most of the structural details are lost in objects this small 
(ISanchez et al.l 120051 ) . For the same reason we also excluded objects touching the edge of 
the map. The perimeter-area log-log plot is shown in Figure H] and the resulting perimeter 
dimension (twice the slope of the best fit) is also listed in Table [TJ The maximum spatial 
scale involved in this calculation can be estimated as \/A max , being A max the area of the 
largest projected cloud in the m ap. We obtained A max = 240 pixel 2 and the pixel size is 6.87 



arcsec (jRosolowsky et al.l 120071 ). then the maximum size is ~ 500 pc. 

Table [TJ summarizes the fractal dimensions calculated for each type of data. The 2D 
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Table 1. Summary of Calculated Fractal Dimensions for M33 



Small spatial scales 3- Large spatial scales 



Sample 




N dat 


n b 


-D/.3D 


-D/,2D 


-D/,3D 


Bright stars 




534 


1.01 ±0.05 


1.0-1.9 


1.93 ±0.03 


2.8-2.9 


Faint stars 




1644 


1.42 ±0.04 


2.2-2.4 


1.89 ±0.02 


2.8-2.9 


HII regions 




617 


1.48 ±0.08 


2.3-2.5 


2.01 ±0.03 


2.9-3.0 


Molecular clouds 


149 






1.98 ±0.04 


2.8-3.0 


CO emission 


map 




1.65 ±0.06 


1.6-1.8 







a For bright and faint stars and HII regions small spatial scale means < 500 
pc and large scale means > 1 kpc. For molecular gas large scale is > 500 pc 
(distribution of clouds) and small scale is < 500 pc (CO map). 

b -D/,2D refers either to the two-dimensional correlation dimension D c (for 
the distribution of stars, HII regions and GMCs) or to the perimeter-area 
based dimension D per (for the CO map). Df^-o is the corresponding three- 
dimensional fractal dimension. 
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Fig. 4. — Perimeter as a function of the area for objects in the CO emission map (Fig. [2]) at 
different intensity levels. The slope of the best linear fit (solid line) is 0.825 ± 0.032. 



fractal dimensions were converted to 3D fractal dimensions using results from previous stud- 
ies. For the distribution of po int objects (stars, HII reg ions and GMCs) we used the results 
for projected disks derived in lSanchez fc Alfarol (I2008, Fig. 1) assuming a flatness (that is, 



a disk thickness to diameter ratio) of 1 (jMa et al.lll998l ) 



For the CO map we used the 
i.e. the maximum distance in 



results obtained in ISanchez et al.l (120051 ) with a "resolution 
the image in pixel units) of N P i X ~ 200 pixels. In any case, the corresponding 3D fractal 
dimension fo r the perimeter-area dimension derived here is very weakly dependent on N pix 
(see Fig. 8 in lSanchez et al.ll2005f ). The procedure of estimating a 3D fractal dimension from 
its corresponding projected value usually increases the associated uncertainties. We have 
included in column labeled -D/,3D hi Tabled] the range of three-dimensional dimensions that 
are compatible with the calculated values of D c or D pcr . It is important to note, however, 
that these relatively large uncertainties do not arise from the method used to calculate Df, 
but from possible random fluctuations in the projection of the original fractals. That is, it is 
not possible in principle (at least using only information from spatial distribution) to disen- 
tangle which of the projected clumps are real and which are the result of chance groupings 
during the projection. 
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Discussion 



The first thing we note is a statistically significant difference between the projected 
dimensions at relatively small spatial scales and those at larger spatial scales. At small scales 
the correlation dimension of the distribution of stars and HII regions is < 1.5 whereas at 
large scales it is > 1.9, and the differences are always larger than the associated uncertainties 
(Table [T|). The spatial scale where this transition takes place is roughly the same for each 
component (~ 500 — 1000 pc). The transition from a smaller correlation dimension to a 



larger one for the young stars in M33 was first reported by lOdekonl (120081 ) . Here we provide 
a more detailed quantification of this transition and observe it for the first time for the 
distribution of HII regions. The transition is not observed for the distribution of GMCs, but 
given the limited number of data points for this sample (N = 149) the lower limit of reliable 
values is higher than for the other objects (r > 500 pc). 

What is the nature of this transition? It is no t a sim ple edge effect. Our algorithm 
corrects for edge effects and, as discussed by lOdekonl ( 120081 ). we would expect an edge effect 
to cause a decrease in the slope on larger scales. On the other hand, some kind of change in 
the measured dimension is expected for projected fractal disks. Given a monofractal of three- 
dimensional dimension -D/,3D, the pr ojected fractal dimensio n -D/,2D of an extracted slice will 
depend on the slice thickness H (see lSanchez fc Arfaroll2008l . for detailed models and discus- 
sion). For very thin slices -D/,2D = -Df,3D — 1 whereas for thick slices -D/,2D = mrn {2, -D/,3d} 
( IFalconerlll990[ ). Then, when calculating the correlation function we should observe a change 
in the behavior around the scale H. For r 3> H the distribution can be seen as a thin slice 
whereas r <C H corresponds to a very thick slice (this latter situation is analogous to a 
simple, direct, projection). However, this is expected to pro duce smaller values of Df 2 r> 
for larger spatial scales (see Fig. 1 in I Sanchez fc Alfarol 120081 ). Thi s kind of transition has 
been observed i n the distribution o f neutral hydrogen in the LMC (jElmegreen et al.l 120011 ) 
and NGC 1058 (IDutta et al.ll2009al ). The two-dimensional power spectra of these galaxies 
flattens at large scales, indicating a smaller fractal dimension; this behavior can be inter- 
preted in terms of a change from 3D fluctuations on small scales to 2D fluctuations on scales 
larger than the disk thickness. In fact, power spectrum analysis h as been suggested as a 



tool to determine t he thickness of nearly face-on galactic gas layer (JElmegreen et al.l 12001 



Padoan et al.l 1200 ll ). This effect contrasts with those presented here for molecular gas and 
young stars, where the correlation dimension is higher on larger scales. 



A different possibility was suggested by iPadoan et al.l ( 1200 ll ) . They argued that there 
must be a physical transition in the statistical properties of the flow close to the disk scale 
height. In turbulent flow the energy is injected at certain spatial scale and then it "cas- 
cades" to smaller scales. But there are many possible energy sources that may be relevant 
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at different levels — for instance, stellar outflows at small scales and galactic shear at large 
scales. A possible consequence may be different distribution patterns at different size ranges. 
Even though the underlying turbulent structure tends to be the same, non-turbulent mo- 
tions acting on galactic scales could modify the final structure at those scales. In other 
words, the power law behavior at small spatial scales would be a direct consequence of the 
self-similar turbulent motions in the medium, but this turbulence is unlikely to extend to 
very large scales, where two-dimensional flows should dominate the dynamics. The transi- 
tion from small-scale three-dimensional turbulence to two-dimensional large-scale motions 
on the disk should occur around the galactic scale height. Interestingly, the behavior we 
observe is that all the fractal dimensions at r > 1 kpc are within a narrow range of values 
(-D/,2D = 1.9 — 2.0, or -D/,3D = 2.8 — 3.0) that are consistent with essentially uniform (ran- 
dom) distributions. Thus, we identify a characteristic spatial scale (around 500 — 1000 pc) 
that separate s two different physica l regim es. This scale is of the order of the typical size 
estimated by lEfremov fc Elmegreenl ( 119981 ) for the larg est star complexes in spiral ga laxies, 
which correlates with the galaxy isophotal radius R25 (jElmegreen fc Elmegreenl I1983I ). The 
transition zone separates regions where coherent star formation is occurring in a turbulent 
medium (stellar complexes) from larger regions that are organized by large-scale galactic 
dynamics. The analysis of the transition in the fractal dimension provides an interesting 
method to determine the typical size of star complexes in a given galaxy. 

Despite the rather large uncertainties for -D/,3D 5 it can be clearly seen that there are two 
separate range of results for the three-dimensional fractal structures at small spatial scales. 
On one hand, bright stars and molecular gas are distributed with -D/,3D ^ 1-9. Young, 
newborn stars should reflect the same conditions of the ISM from which they were formed, 
which is roughly consistent with this result. On the other hand, faint stars and HII regions 
have -Dj j3 d ~ 2.2 — 2.5. These higher fractal dimension values can be in terpreted in terms o f 
evolutiona ry effects. The clu mpy distributions of young stars in LMC ( JBastian et al.l 120091 ) 
and SMC (IGieles et al.ll2008l ) evolve towards smoother d istributions as stellar ages increase. 
The same is true for the stars clusters in these galaxies ( IBonatto fc Bicall2010f ). It has been 
shown that the brightest HII regions in spiral galaxies (which reflect, in a first approximation, 
the initial distribution of star-fo rming sites) tend to be d istributed in more clumpy patterns 
than the low-brightness regions ( ISanchez &: Alfarol 120081 ). The higher fractal dimensions for 
faint stars and HII regions at small scales are probably due to this kind of evolutionary effect. 

We can compare the small-scale three-dimensional structure of the interstellar medium 
in M33 with that of our own galaxy. The spectral index 7 of the power spect r um in the Milky 



Way lie s in the range 2.8 < 7 < 3.2 in 2D maps (JElmegreen fc Scald 120041 ) . iBrunt fc Heyer 



( I2002bl ) systematically measured the energy spectrum in several CO emission maps of the 
outer Galaxy and obtained /3 = 2.17 ± 0.31, significantly higher than the Kolmogorov spec- 
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trum P = 5/3 ~ 1.67. The corresponding dimension of the projected contours is D pcr ~ 1.4, 
very similar to the avera ge value found from direct measurements of D per in the Galaxy 
(jElmegreen &: Scaldl2004l ). In a previous work ( jSanchez et al.ll2007bl ) we used various molec- 
ular emission maps to study the fractal properties of Ophiuchus, Perseus, and Orion, and 
found that the dimension is always in the range D pei ~ 1.30 — 1.35. See Appendix A for a 
discussion of the definitions and relationships among 7, (3, and D 



per- 



The 3D fractal dimension of the distribution of molecular gas for M33 (-D/,3D — 
1.8, Table [H), meanwh i le, is m uch smaller than for the Milky Way (-D/,3D — 2.6 



1.6- 

- 2.8, 



sec 
to /3 



Sanchez et al.l 120051 l2007bl ). The CO emission map va lue of D 



per 



1 . 65 corresponds 



1.7, very close to the Kolmogorov value of 5/3q Elmegreen et al.l ( 20031 ) found, 
from azimuthal scans of M33, that the power spectrum of optical emission has an index 
even smaller than the Kolmogorov value, although the uncertainties are relatively high. 
That is, the boundaries of the projected perimeters in M33 are on average more irregular 
(higher D per values) than those of the Milky Way because molecular clouds in M33 exhibit a 
more fragmented structure (smaller -D/,3D values). The fact that -D/,3D ^$1-9 for the three- 
dimensional distribution of very young s tars is consistent with the small fractal dimension of 
the CO gas. In addition. iBastian et al.l ( 120071 ). using the minimum spanning tree technique, 
found that the number of objects increases with the radius as iV ~ R a with a = 1.63 — 1.96. 
Random sampling in perfect hierarchical fractals yields N ~ R D f so that, in general, their 
results are consistent with the relatively small fractal dimensions we are reporting here. 
Whether there are (or not ) othe r differences with GMCs in our Galaxy remain an open 
problem. lEngargiola et al.l (120031 ) found that the mass spectrum of GMCs in M33 follows 
a power-law distribution (according to a hierarchical picture) but with a slope considerably 
steeper than that found in the Milky Way. However, ISheth et al.l (120081 ) concluded that if 
the same cloud identification algorithm and analysis technique are applied, and if resolution- 
dependent effects are taken into accou nt, then there seem to be no differences in the GMCs 
properties between these two galaxies. iBolatto et al.l (120081 ) analyzed in a consistent manner 
GMC properties (sizes, velocity dispersions, and luminosities) in several galaxies (including 
M33) subject to a wide variety of physical conditions and found only small differences with 
the properties of GMCs in the Milky Way. However, it is important to keep in mind the 
difficulty in comparing observations of an external galaxy with observations of the Milky 
Way. Molecular clouds from large-scale surveys in the Milky Way can be blended in space 
and velocity, making it troublesome to extract information on their internal structures. 
Blending effects are reduced in external galaxies such as M33 that are seen nearly face-on. 



3 This is only a reference value. We are not saying that turbulence in M33 actually behaves according to 
Kolmogorov's description. 
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Moreover, the internal structure of GMCs in the Milky Way is measured for clouds in the 
solar neighborhood, at an average Galactocentric distance of about 8 kpc. Instead, the 
average internal structure of clouds in M33 was calculated in this work for the molecular 
map of the central region, at galactocentric distance less than ~ 2 — 3 kpc. Obviously, more 
detailed studies are needed to clarify this point but, in this work, we have found a significant 
and remarkable difference in the internal structure of GMCs between M33 and the Milky 
Way. 



5. Conclusions 

In this work we measure the clustering in the distribution of stars, HII regions, molecular 
gas, and individual giant molecular clouds in M33 as a function of spatial scale. We identify a 
transition in the clustering strength from a very clumpy (fractal) structure at relatively small 
scales to almost uniform patterns at larger scales. The spatial scale for this transition is in the 
range ~ 500 — 1000 pc, independently from the type of object considered. We argue that this 
characteristic scale separates two physical regimes, one where small-scale turbulent motions 
generate self-similar structures and another dominated by large-scale galactic dynamics. 

The existence of this transition implies that care must taken in calculating a single 
value for the fractal dimension over a large range of scales. One must take into account not 
only edge effects on large scales and sampling/resolution effects on small scales, but also 
the possibility of physical changes on intermediate scales. Ideally, this could be done using 
the differential form of the correlation function combined with an edge-correction algorithm; 
this way, changes with scale are not carried into the calculation of the dimension for larger 
scales. However, for small data sets (even the larger data sets presented here, with on the 
order of 1000 objects), the differential form is rather noisy. 

At small spatial scales, bright stars and molecular gas in M33 are distributed with nearly 
the same three-dimensional fractal dimension. This result is consistent with the idea that 
newborn stars follow the same patterns of the ISM from which they were formed. Faint stars 
and HII regions exhibit higher fractal dimensions possibly as a consequence of evolutionary 
relaxation. 

Interestingly, the three-dimensional fractal dimension of molecular gas in M33 (-D/,3D ^ 
1.9) is significantly smaller than in the Milky Way (-D/,3D — 2.7 ± 0.1). This means that the 
interstellar medium in M33 is on average more fragmented and irregular than in our own 
galaxy. This result may have important implications for the study of the physical processes 
that determine the ISM structure. 
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A. Appendix 

Interstellar turbulence can be characterized in several different ways. Two commonly 
used measures are the energy spectrum E(k) ~ \k\~@ and the power spectrum P(k) ~ |&|~ 7 
for wavenumber k. The energy spectrum is the integral over all directions of the power 
spectrum, E{k)dk = J P(k)dfl, so that with this nomenclature the scaling exponents /3 and 
7 are related through fjBrunt fc Heyerll2002al ) 



where E is the euclidean dimension of the image. For example, for a dissipationless cascade 
of energy through an incompressible fluid the Kolmogorov energy spectrum is /3 = 5/3 
and the power spectrum index is 7 = 11/3, 8/3, or 5/3 in 3D, 2D or ID, respectively 
felmegreen fc Scald I2004} ). 



The structure of molecular clouds can also be described as a fractional Brownian motion 



(fBm) structure (IStutzki et al.lll998t iMiville-Deschenes et al.l 120031 ). The properties of fBm 
clouds depend on a single parameter, the so-called Hurst exponent H. This parameter 
character izes the self-similar structure a nd it is related to the corresponding power spectrum 
through (JMiville-Deschenes et al.ll2003l ) 



~f = 2H + E . 

The iso - intensity contours of the i?-dimensional image of a fBm-cloud is given by D lso = 
E — H (IStutzki et al.lll998l ). which corresponds to the perimeter- area based dimension D per 
for two-dimensional maps. In this last case (E = 2) we can write the equations relating 
these commonly used parameters as 



D 



per 



H 
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D per = (5 - 0)/2 , 

and 

D pcr = (6 - 7 )/2 . 

For very rough structures (H = 0), we have (5 = 1, 7 = 2 and _D pcr = 2, while for very 
smooth structures (H — 1) we get /3 = 3, 7 = 4 and -D per = 1. 
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